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Abstract 

We compute the asymptotic determinant of the discrete Laplacian on a simply-connected rectilinear 
region in R 2 . Specifically, for each e > let H e be the subgraph of eZ 2 whose vertices lie in a fixed 
rectilinear polygon U. Let N(H^) denote the number of vertices of H t and B(H e ) the number of vertices 
on the boundary (the outer face). Then the log of the determinant of the Laplacian on H e has the 
following asymptotic expansion in e: 

^N(H t ) + l ^^B m - J, 2 (e, U) + 0(1) 

where G is Catalan's constant and r2(e, U), which is 0(log i), is the Dirichlet energy of a certain canonical 
harmonic function h on U. 

As an application of this result, we prove that the growth exponent of the loop-erased random walk 
in Z 2 is 5/4. 

Resume 



Nous calculons le developpement asymptotique du determinant du Laplacien discret sur une region 

le application, nous mc 
aleatoire laplacienne dans Z 2 est a croissance n 5 ^ 4 . 



rectilineaire de R 2 . Comme application, nous montrons que l'esperance de la longueur de la marche 



1 Introduction 

The determinant of the Laplacian on a graph arises in two related statistical mechanical models, the 
uniform sp ann ing tree model and the two-dimensional lattice dimer model. For the spanning tree model, 



Kirchhoff [ Kii ] is attributed with showing that the product of the non-zero eigenvalues of the Laplacian 



on a finite graph is the sa me a s the number of spanning tree s on th at graph. For the two-dimensional 



dimer model, Temperley [rem], based on work of Kasteleyn [Kasl , showed that the number of dimer 
coverings of certain subgraphs of 1? can be computed by the determinant of the Laplacian on related 
graphs. Precise estimates on these determinants provide important information about these mode ls, in 



particular allowing one to compute certain critical exponents and correlation functions [DD, Kenl|. 

In this paper we compute the asymptotic expansion of the determinant of the Laplacian on a special 
family of graphs: subgraphs of Z 2 which are approximating rectilinear polygons (a polygon is rectilinear 
if its sides are parallel to the axes). Our main motivation is not to study the Laplacian in itself but rather 
to study both the dimer (domino tiling) model and the uniform spanning tree model. For this reason we 
use the language of domino tilings. (Domino tilings are tilings with 1x2 and 2x1 rectangles.) 
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Temperley [ Terr ] gave a bijection between the number of spanning trees of a subgraph H of Z and 
domino tilings of a polyomino P = P(H) constructed from the superposition of H and its dual. We call 
a polyomino Temperleyan if it arises from a graph H by Temperley's construction. Such polyominos 
have a simple description: let Q be a polyomino such that each side between a concave and convex 
corner has even length, and each side between two concave or two convex corners has odd length. Such 
a polyomino has odd area; let P be obtained from Q by removing one lattice square at some convex 
corner. Then P is Temperleyan. 

Theorem 1 Let U C R 2 be a rectilinear polygon with V vertices. For each e > 0, let P e be a Temperleyan 
polyomino in eZ 2 approximating U in the natural sense (the corners of P t are converging to the corners 
of U). Let A t be the area and Perim t be the perimeter of P t . Then the log of the number of domino 
of P € is 

coA t c\Perim e n ( 1 \ 

+ c 2 (e) log - + c 3 (U) ) + c 4 + o(l), (1) 



e^ e 48 \ e 

where co = ^-, G=l— + — ... is Catalan's constant, ci — ^ + los( -^ 2 ~ 1 - > , C4 is a constant 
independent of U and 02(e) log f + cz(U) is the e-normalized Dirichlet energy of the limiting average 
height function on U (see definitions below). The term — jgC2(e) is of the form 

1 V — 4 

(1 + ERR(e)) 
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where ERR(e) is o(l). 

Corollary 2 Under the above hypotheses on U, for each e > let H € be the subgraph of eZ 2 whose 
vertices are in U. Let N(H e ) be the number of vertices in H t and B(H t ) the number of edges of el? on 
the boundary of H t . Then the log of the determinant of the Laplacian on H e is 

fN(H e ) + l °^-^ B(H t ) - I ( C2(e) log \ + c 3 (U)) + cs + o(D, 

where C2,C3 are as in Theorem^ and C5 is another constant independent of U . 

The limiting average height function has the following description. Let 60 £ dU be a base point. For 
x € dU define uo(x) to be the total turning (in radians) of the boundary tangent on the boundary path 
counterclockwise from feo to x (the function «o has jump discontinuities at each corner of U and at feo)- 
The limiting average height function is the harmonic function on U whose boundary values are — uq. The 
e-normalized Dirichlet energy is by definition the Dirichlet energy contained in the complement of the 
e-neighborhoods of the jump discontinuities. 
Remarks 

1. In the case U is a r ectang le, formula (|l|) follows from the formula for the exact number of tilings 



computed by Kasteleyn [Kasl and Temp erley and Fisher |TF| (the asymptotic expansion of which was 
computed by Duplantier and David ppj); see Proposition |16] below. 

2. The leading term in th e above formula, involving the constant Co, essentially follows from work of 
Burton and Pemantle ]BP| : They constructed a measure /x of entropy Co on the space X of domino tilings 
of the plane and proved that it was the unique translation-invariant measure of maximal entropy on X. 
Furthermore they proved that for regions of the type used in the theorem, the entropy (the coefficient of 
e" 2 in (0)) is c Q A. 

3. Our boundary conditions give rise to a correction to the number of til ings w hich is only exponential 



in the length of the boundary (the 'perimeter' term in the theorem) . In | CKF ] , on the contrary, it was 



shown that in some sense "most" other boundary conditions have a larger effect, giving a smaller entropy 
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Co, and making the local densities of configurations vary throughout the region. So both the 'Area' and 
'Perimeter' terms in the theorem depend strongly on our choice of boundary conditions. 

4. Note that if two regions have the same area, perimeter and number of vertices then the log of the 
number of domino tilings differs by a constant in the limit, that is, the ratio of the number of tilings is 
tending to a constant as e — > 0. This constant depends on the shape of the regions and can in principle 
be computed explicitly. 

5. There has been much work done on the 're gular i zed' deter minant of the continuous Laplacian, and 
the asymptotic distribution of its eigenvalues |Kac, VIS, OPS]. We have not attempted here to make 



any c o nnec tion between these two subjects, although there is a lot of evidence for a connection, see e.g. 



OPS 
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6. As noted above, Temperley |Tem[ gave a bijection between the set of spanning trees of a subgraph 
of Z 2 and the set of domino tilings of a related polyomino. The corollary follows from the theorem by 
applying this bijection (see section |i| for the definition) : the graph H t gives rise to a polyomino P e of 
area A t = e 2 (AN(H e ) — B(H e ) — 4) and perimeter Perim e = e(2B(H,,) + 4). Plugging these values into 
(Q) gives the formula in the corollary. 

7. The function ER R(e) i s unknown although it seems possible that it could be computed using Toeplitz 
determinants, as in [MW . 



Part of the motivation for proving Theorem [l] is to validate a certain heuristic, which attempts 
to explain how the presence of the boundary affects the long-range structure of a random tiling. In 
particular it at tempt s to explain how the boundary affects the densities of local configurations far from 
the boundary [DMB|. We call this heuristic the 'phason strain' principle. 

The h euris tic is as follows: the boundary causes the average height function of a tiling (see definition 
in section 2.3) to deviate slightly from its entropy-maximizing value of 0. At a point in the region where 
the average height function has nonzero slope, the "local" entropy there is smaller than the maximal 
possible entropy, by an amount proportional to the square of the gradient of the average height function. 
The system behaves in such a way as to maximize the total entropy subject to the given boundary values 
of the height function, and the resulting average height function is the function which minimizes the 
(integral of) the square of its gradient. That is, the average height function is harmonic. This "explains" 
the terms 02(e) log - + cs(U) in Theorem [lj 

Unfortunately the const ant -p appearing in Theorem [l] is different from the expected value of the 
local entropy as deriv ed in ICKP | (where a rigorous version of the phason strain principle is proved in a 
different context): in CKP| the entropy as a function of slope is shown to have the expansion 



ent{s, i) 



ent(O.O)- -^(s 2 +t 2 ) 



+ 0(terms of order > 3), 



and where (s, t) are the p artial derivatives of the height function. We may conclude from this discrepancy 
that the computation in [CKP| can not be re fined t o obtain asymptotics of the same precision as Theorem 
^ above (when applied to the present case, [3KP only gives the leading term in (Q)). In particular the 
phason strain principle can not be considered valid in this context. 

The techniques used to prove Theorem |l| can be applied to the uniform spanning tree model as well. 
Indeed, as mentioned abo ve, there i s a clos e connection between the spanning tree process on 1? and 
the domino tiling model Tern, BP, KPW[ . Many properties of spanning trees on 1? translate into 



computable properties of dominos. We study here one particular property of a uniform spanning tree: 
the distribution of the (unique) arc between two fixed points. The relevant question about tilings is to 
count the number of tilings of a region with a hole (single square removed). To estimate this number, we 
use the technique of Theorem we cut the region apart up to the hole, and then sew it up again in such 
a way as to remove the hole. In this way we prove the well-known conjecture that the expected number 
of points on the tree branch within distance N from the origin grows like N 5/4 . In fact we prove more: 
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Theorem 3 On the uniform spanning tree process on N x Z, the expected number of vertices on the 
branch from (0, 0) to oo which lie within distance N of the origin is 
of a vertex (x, y) = re' e to be on the branch from (0, 0) to oo is 

r -3/4(l + /M ) cosW l/4 (1 + o(1)) 

where f(r) is o(l) as r — > oo. 

The bran ch in a uniform spanning tree has the same distribution as the loop-erased random walk (LERW) , 



see [Pern | . So this proves that the growth exponent of the loop-erased random walk is |. This value 



of the exponent d has been conjectured by physicists for some time [GB, Maj], using arguments based 



on conf orma l field theory and the assumption of conformal invariance of the "scaling limit" of the walk. 
Lawler |Law[| had previously given the bounds 1 < d < |. 

Here is an outline of the paper. Section |^ gives the definitions and background. Most of the back- 
ground comes from [Kenl| and [ Ken2J : loc al properties of dominos can be found in [Kenl|, and the 
conformalproperties can be found in |Kcn2| . In Section ^, we state two lemmas and use them to prove 
Theorem hi Specifically, the theorem is proved by cutting up a rectilinear polygon into rectangles and 
using the known formula for the number of tilings of a rectangle. Lemma ^| determines how the number 
of tilings changes as you are making a single cut, and Lemma [?] relates this change to the change in 
Dirichlet energy of the average height function. The next two sections are devoted to the proofs of the 
lemmas. Section |^ recalls the formula for the number of tilings of a rectangle and proves the formula 
of Theorem [l] in this special case. Section Q discusses the connection of domino tilings to loop-erased 
random walk and proves Theorem [| 

We kindly acknowledge Oded Schramm, Wendelin Werner and Bertrand Duplantier for helpful dis- 
cussions, and thank Oded Schramm and Russell Lyons for proofreading. 



2 Definitions and Background 

2.1 Temperleyan polyominoes 

By the grid Q we mean the graph whose vertices are 21? and edges join all pairs of vertices at distance 
2. A lattice square is a face of Q. A simply-connected subgraph of Q is a set of vertices and edges 
of the grid which is the 1-skeleton of a simply connected union of (closures of) lattice squares. 

Let H be a finite simply-connected subgraph of Q and b £ H a fixed vertex adjacent to the outer 
face of H. We associate to H a new graph P' — P'(H) = P'(H,b) as follows. There is a vertex of P' 
for each vertex, edge and face of H, except for the outer face and the vertex b. Two vertices u\,u% of 
P' are connected by an edge in two cases: m , ui come from an edge e and a vertex v of H (and v is on 
the edge e), or ui,U2 come from a face / and edge e of H (and e is part of the boundary of /). In other 
words, P'(H) is the "superposition" of H and its planar dual H' , except that we discard vertex b of H 
and the outer vertex of H' . Let P{H) denote the polyomino in whose dual (not including the vertex 
for the outer face) is P'(H). See Figure [ljfor an example. 

Except for the outer face, the faces of P axe squares, and come in four types, Bo, Bi, Wo, Wi. The 
squares in Bo are those coming from vertices of H; the squares in B\ are those coming from faces of H. 
The squares in Wo (respectively W\) are those coming from horizontal (resp. vertical) edges of H. The 
'B' and 'W stand for 'black' and 'white' coming from the checkerboard coloring of P. See Figure |l| We 
assign colors B p,B\ , Wo, Wi to vertices of P' corresponding to the colors of the faces of P. 



Temperley |Tem] showed that there is a bijection between spanning trees of H and perfect matchings 
of P'(H) (a perfect matching, or dimer covering, is a set of edges such that every vertex is contained 
in a unique edge). Perfect matchings of P' correspond to domino tilings (tilings with 1x2 and 2x1 
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Figure 1: The graph H and polyomino P(H). The squares of Bo are black, those of B\ are in grey. The 
vertex b of H is in grey, and corresponds to the missing corner of P. 



rectangles) of P(H). A polyomino (a union of lattice squares bounded by a simple closed curve) is said 
to be Temperleyan if it is of the form P(H) for some simply-connected subgraph H of Q. Note that if 
H is a simply-connected subgraph of eQ for some e > 0, then P(H) will be a polyomino in 

The lattice square missing in P which corresponds to point b £ H is called the base square of the 
polyomino. For the graph P' dual to P, b is the base vertex (it is not a vertex of P'). By Temperley's 
bijection, the number of tilings of P is independent of choice of vertex b £ H as long as it is on the 
boundary. If P is a Temperleyan polyomino then by H — H(P) we mean the unique associated simply 
connected subgraph of Q for which P = P(H). 

These definitions also apply to infinite graphs H, the only difference being that in this case we may if 
we like choose b = oo, which means we do not remove any lattice square from P(H). All infinite graphs 
we deal with in the sequel have the property that near to infinity the boundaries are straight, that is, 
the boundaries have no corners outside some fixed large radius. This is to avoid certain convergence 
problems later. Two important examples of infinite Temperleyan polyominos are the whole plane P(l?) 
and the half plane P(N x Z). Both of these have b — oo. 



2.2 Conformal properties 

The results of [ Ken2 apply to Temperleyan polyominoes. Let P be a Temperleyan polyomino with dual 
graph P'. We assign weights to the edges of P' so that a horizontal edge has weight I if its left vertex is 
white, —1 if its left vertex is black; a vertical edge is weighted i — v^— T if its lower vertex is white, — i if 
its lower vertex is black. Thus the weights around a white vertex are 1, i, — I, —i in counterclockwise order 
starting from the edge leading right; around a black vertex these weights are — 1, — i, The adjacency 
matrix of the graph P' with these weights is call ed th e Kasteleyn matrix Kp of P. Its determinant is 
the square of the number of domino tilings of P |Kas2|. The inverse of the Kasteleyn matrix is called the 
coupling function Cp(-,-) of P. The probability of a configuration of dominos occurring in a random 
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tiling is the absolute value of the determinant of a submatrix of the coupling function matrix [Kenl . 
The coupling function has a number of important properties which we list here. 



2.2.1 Combinatorial properties of the coupling function 

First, for vi,V2 two vertices of P , we have Cp(vi,V2) = Cp(v2,vi) and if vi and V2 are both black or 
both white then Cp(v 1,^2) = 0. Therefore we will always take the first variable of Cp to be a white 
vertex and the second variable to be a black vertex. 

The coupling function has a concise description in terms of the Green's function on H . Let G be 
the Green's function on the graph H, that is, G is the function on H x H which satisfies AG(x, y) = 
S x (y) — Sb(y), where A is the Laplacian with respect to the second variable (and recall that b is the base 
vertex). Here A/(u) = 4/(f) — f(v + 2) — f(v — 2) — f(v + 2i) — f(v — 2i), except at a boundary vertex, 
where A/(f) is the degree of v times f(v) minus the sum of the neighboring values. (Here v ± 2, v ± 2i 
refer to the four neighbors of v in the graph H. Note that these vertices are at distance 2 in P' .) As 
stated the function G is only well-defined up to an additive constant. We fix the constant by setting the 
function to be zero when y — b. For x,x' any two vertices of H, the function L(y) := G(x,y) — G(x',y) 
satisfies AL(y) = 5 x (y) — 8 x i(y) and L(b) — 0. 

Let x and x' be adjacent vertices of H; let /, /' be the faces of H adjacent to the edge xx' , with /' 
on the left as the edge is traversed from x to x' . Let L be the function on the faces of H which is the 
harmonic conjugate of L(y) in the sense that for any edge e = v\V2 of H which is not the edge xx' we 
have L{v2) — L(vi) = L(fi) — L(f2), where /1 is the face to the left of the edge e (when e is traversed from 
vi to V2) and /2 is the face to the right. If we define L to be zero on the outer face then it is uniquely 
defined and harmonic except at / and /'. Moreover AL(-) = Sf(-) — 5f/(-), where A is the Laplacian 
on the dual H' of H . The function L(z) can also be written G(f, z) — G(f' , z), where G is defined by 
AG(f,z) — Sf(z) — S (z) and G(f,o) = 0, o referring to the outer face. 

We now have the following description of the coupling function in terms of these Green's functions. 
Suppose vi £ Wo- Then 

r , s_( G(vi + 1,112) - G(vi - I.V2) if V2 G Bo 

Op(V!,V2)-^ 1 (g(„ 1+1j1;2 )_g(„ 1 _ 1j „ 2 )) i £ Va€Bl . ( ' 

Note that when V2 G -Bo, Cp(v\, V2) = G(vi + l,Ua) — G(vi — l,i>2) makes sense since both v± ± 1 and 
V2 axe in Bo (hence vertices of H). When V2 £ Bi rather, Cp(vi,V2) = i(G(v\ + i,V2) — G(v\ — i,V2), 
makes sense since vi ± i and V2 are in B\ (faces of H). Similarly, if v\ £ Wi we have 

r / \ = / G(vi + l,v 2 ) - G(vi - l,v 2 ) if V2 G B± , . 

(Vl ' 2> ~\ -i(G(vi+i,Va)-G(vi-i,V2)) ifv 2 £B . 1 ' 

This description of the coupling function follows from the fact that K*K is the Laplacian on H , and 



also acts as the Laplacian on the dual of H, see [Ken2 



An important case of the above formula for Cp is when vi corresponds to an edge on the boundary 
of H , so that one of V\ ± 1 or Vi ± i corresponds to the outer face of H. Suppose for example that v\ + 1 
is the outer face of P'; then 5 vi +i(v2) = by definition (for V2 a vertex of B\, that is, a face of H which 
is not the outer face) and so the term G(v\ + 1, V2) can be ignored in the above formula. 



2.2.2 Asymptotic properties 

Let U be a rectilinear polygon in C. Fix a base point 60 G dU. Let {P e } e> o be a sequence of Temperleyan 
polyominoes P t C eZ 2 , approximating U as e — > in the following sense. The P e are rectilinear with the 
same number of corners as U, one corner converging to each corner of U. Furthermore the base points 
b e G P t converge to 60. 
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In [ Ken2| we proved the following. Let C e = Cp t . Under the above convergence hypotheses, the 
rescaled coupling functions \Ct converge to a pair of complex- valued functions Fo(v, z) and Fi(v,z) 
which are meromorphic in z, in the following sense. If {u e }, {v t }, {w £ }, {i e } with u e ,v e ,w e ,x e G P e are 
four sequences of vertices of type Wo, Wi, Bo, Bi respectively, converging to respectively u,v,w,x G U, 
then 

lim — C e (u e ,w e ) = Re_Fb(it, w) 

e— »o £ 

lim —C e (u e , Xe) = ilnxFnfw, x) 

lim -CAvt.wA = ReFi(v,w) 

lim —Ce(v e ,x e ) = Hm.Fi(v, x). 
e->o e 

The functions Fq and Fi are defined by the following properties. 



Proposition 4 [ Ken2 , Theorem 13] For each fixed v £ U the function Fo(v, z) has a following prop- 
erties: 

1. it is meromorphic as a function of z G U 

2. its imaginary part vanishes for z G dU , 

3. it has a zero at z = bo, 

4- it has a a simple pole of residue — at z — v, and no other poles on U . 
Similarly, for each fixed v G U the function Fi has the following properties: 

1. it is meromorphic as a function of z G U; 

2. its real part vanishes for z G dU, 

3. it has a zero at z — bo, 

4- it has a simple pole of residue — at z — v, and no other poles on U. 
Furthermore, Fo and Fi are the unique functions with these properties. 

Define the functions F+ = F+ := Fo + Fi and F- = FE := Fq — Fi. These functions are easier to 
work with since they transform as nicely under conformal mappings: F+(v, z)dv is a meromorphic 1-form 
and F-(v, z)dv is an antimeromorphic 1-form. 



Proposition 5 [Ken2, Proposition 15] The function F+(v,z) is analytic in both variables. The 
function F-(v,z) is analytic in z and anti-analytic in v. IfV is another marked region and f: V — > U a 
conformal isomorphism sending the base point ofV to the base point ofU, then we have the transformation 
rules 

F^(v,z) = f'(v)F^(f(v),f(z)) (4) 
FY(v,z)=T^)FY(f(v),f(z)). (5) 

When U = C with bo = oo we have F + (v,z) = n ^- v ) ant ^ ^-(v,z) = 0. When U is the right 
half-plane RHP — {z\Ke(z) > 0} with bo — oo we have 

F+(V,Z) = — r-, F-(V,Z) = j^- r. (6) 

This and the above transformation rules determine F+ and F- (and hence Fo,Fi) on any simply- 
connected region. 
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2.3 Average height function 



Recall Thu] that the height function of a domino tiling is an integer- valued function on the vertices of 
the dominos; it is well-defined up to an additive constant. After fixing its value at some vertex, it is 
defined by the property that on an edge viv? which is not crossed by a domino, the height difference 
h(v2) — h(vi) is +1 if the square to the left of the edge V1V2 (we mean, the square containing edge V1V2 
and on the left when traversing V1V2 from vi to V2) is black, and —1 if this square is white. See Figure 
^| Note that the height function along the boundary is independent of the tiling. 
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Figure 2: Some of the heights in a domino tiling. 



Let U be a rectilinear polygon with base point 60 and P e a sequence of Temperleyan polyominos 
converging to U as in the previous section. Suppose for simplicity that all the P t contain a fixed vertex 
vo, say at a corner of P t . Let h e be the average height function of P t , that is, for any vertex v £ P e , 
h t (v) is the average height of v over all domino tilings of P e , where the height at vo is taken to be zero. 
The limiting average height function h of a random domino tiling of U is by definition the limit as 
e — > of the functions h e : take x £ U, with x dU and let x t £ P t converge to x; then 

h(x) := lim h t (x t ). 

For x g dU b ut not at a corner, h(x) is defined by continuity from values of h in the interior. 



In [Kcn2] we showed that this limit exists and has a simple expression in terms of the function _F+. 
Let F+{v,z) = F+(v,z) — ant ^ ^ F+( v ) = li m z-»u F+( v i z )- Then the limiting average height 

function h on U is given by the complex line integral 

h(v) = 21m / F'+(u)du. 

JVQ 

The choice of zero no is immaterial since the height is only defined up to an additive constant anyway. 
Note that from (^|), on C or on the right half-plane we have F+(v, z) = n ^_ v ^ so F+ = and therefore 
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h(v) is constant, as expected. As another example, the map z i±£ maps the unit disk to the 
right half-plane, mapping 1 to oo; therefore using (^), on the unit disk with basepoint z = 1 we have 
F+(v,z) = n 2(1 ~ z) - and Fl(z) = -r^, so that h(v) = --Im P 

' / 7r(l — Uj(z — 1>) + V / 7r(l — Z) ' V / 7T J tip 1 _ W 

Note that in general h(v) is harmonic, being the imaginary part of an analytic function. 

For v on the outer boundary of U , h(v) — h(vo) is the total turning (in radians times 2/tt) of the 
boundary tangent on the path counterclockwise (cclw) around the boundary from Vo to v. However when 
passing the basepoint bo the height drops by 4. Thus a full cclw turn contributes +4 — 4 = to h. The 
function h is constant on the straight edges, and there is a discontinuity at each corner of U, where h 
changes by ±1 according to whether the corner is a left turn or right turn (if the basepoint bo is at a 
corner then the change would be —3 or —5 accordingly). These boundary values can be understood in 
a sense using Figure the polyomino in this figure can be thought of as approximating a rectilinear 
octagon U , where the height on the lower boundary of U is —1/2 (the average of —1 and 0, the alternating 
heights on P e ), the height on the right-most boundary is 1/2 (the average of and 1) and so on. Since a 
harmonic function is determined by its boundary values this provides a simple description of h in terms 
of the turning of the boundary tangent. 

2.4 Dirichlet energy 

Recall that the Dirichlet energy of a harmonic function ft on a region U is given by 

E(h)= [[ \Vh\ 2 dxdy = <f hdg 

where g is a harmonic conjugate of h, that is, g is a harmonic function so that h + ig is locally an analytic 
function of x + iy. 

We will be interested in the Dirichlet energy of harmonic functions which are the limiting average 
height functions on rectilinear polygons U. In particular their boundary values have a finite number of 
jump discontinuities (at the corners); unfortunately in such a case the Dirichlet energy is infinite. To 
avoid this difficulty, for each sufficiently small <5 > we define the 5-normalized Dirichlet energy 
Es(h) as follows. Remove a S- neighborhood of each x G dU for which the harmonic function h has a jump 
discontinuity. Let U' be the region U without these neighborhoods. The (5-normalized energy Es(h) is 
simply the integral of |V/i| 2 over U' . 

If U is unbounded and if h has a jump discontinuity at oo, we remove the neighborhood of oo consisting 
of points | z | > j, and compute the energy on the remaining region as before. 

Here we will illustrate with an example. Let h be the bounded harmonic function on the upper half 
plane which has value on the x-axis to the right of the origin and 1 on the a;-axis to the left of the origin. 
Then h(z) = -Imlog(z). The harmonic conjugate to h is g(z) = — -Relog(z). The normalized Dirichlet 
energy is the integral over dU' of hdg, which can be broken into four parts: the integral from — S^ 1 to 
— <5, the integral around the half circle of radius 5, the integral from 8 to <5 _1 , and the integral around 
the half-circle of radius <5 _1 . The third of these integrals is zero since h is zero on the positive a;-axis. 
The second and fourth are zero since dg is zero on circles about the origin. So the only contribution is 
from the first integral which gives 

E s (h)= [ 6 ldg = g(-d) - g(-l/S) = - log i 
J-i/s k 5 

In a more general situation dg will not vanish along the boundaries of the ^-neighborhoods of the 
discontinuities, but dg will still be O(S) there, as we now show: Let h be the average height function on 
a rectilinear polygon U. Take a convex corner of U, translate and rotate U so that the corner is at the 
origin and is bounded by the positive axes. Without loss of generality (after a linear scale) we suppose 
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that near the origin ft is on the a>axis and —1 on the y-axis. Then ft is the imaginary part of an analytic 
function h(z) on U whose expansion at is of the form 

2 2 
— logfaia + a.2z +...) = log(z) + /?() + fi\z + fez + . . . . 

7T 7T 



In particular when z = Se 10 , we have 

dg{z) 



= -Re(/3i5ie lS + 0(S 2 )) = 0(S). 



d9 7T 

A similar argument works at a non-convex corner. 

Since ft has a standard form near each of its jump singularities, the ^-normalized Dirichlet energy has 
a very simple dependence on 8. Recall that at a convex corner the height function changes by +1 when 
moving cclw around the boundary. So the height function is of the form of that of the above example. 
If we change 8 to a smaller 8' , the change in energy at that corner is 



[ lS -dg = --\og(8') + -log 5 + 0(8). 

JiS T 77 



Thus the dependence on 8 at a convex corner is — log 4 + 0(8). 

Similarly we can do the calculation at a concave corner. Move and rotate the corner so that it is 
bounded by the positive x-axis and the negative y-axis. Up to an additive constant ft. is on the x- 
axis near the origin and 1 on the negative y-axis. So ft is the imaginary part of an analytic function 
h(z) = log z + cto + aiZ + oi2Z 2 + . . . . Now if S changes to the smaller 3' , the change in energy is 

-id' 



[ 1 d g = ~\oz{5') + l-\og6 + 0{5). 
J-iS 3tt 3tt 



Thus the dependence on 8 at a concave corner is ^ log j + O(S). 

If a corner contains the base point, then if it is a convex corner the height changes by —3 rather than 
+1. the dependence on 8 is therefore 9 times that for a normal convex corner, or — log 4. So the fact 
that the corner contains the base point adds — log | to its "local energy" . If the corner is concave the 
height changes by 5 rather than 1, and so the local energy is |^ log | rather than j- log j, which is also 
an addition of ij- log 4 . If the base point occurs along an edge, the height change is 4 and the energy 
associated is again i£ log i . 

Now from these calculations, for any rectilinear region U we can immediately compute the dependence 
of the 5-energy on 8, in terms of the number of vertices. A (simply-connected) region with V vertices 
has (V — 4)/2 concave vertices and (V + 4)/2 convex vertices, one of which we may take to be the base 
vertex. So the dependence on 8 of its 5-energy is 

( 2 fV-4\ 2fV + 4:\ 16V 1 
^ + {^{—) + ^{—) + -) l0 Z-8 = 

f4(V-4) 24\, 1 

; + — log 7 + 0{8). 



1^3^ ty J 8 

Note that —tt/48 times the above <5-energy gives the logarithmic terms in Theorem [j] (if we replace 8 



with e) 
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3 Cutting a lattice region into rectangles 

Here we prove Theorem jjj. Let U be a rectilinear polygon with basepoint bo £ dU. Let 70 = 70 (t) be 
a straight (horizontal or vertical) unit speed path in U from dU to dU, which avoids 60 and does not 
touch dU except at its endpoints. 



Let P e C el? be a Temperleyan polyomino approximating U as described in section 2.2.2 Let 7 e be 
a strip of width e of lattice squares of P e , lying within 0(e) of 70 and traversing P e from the boundary 
to the boundary, avoiding the base square of P e . Furthermore we require that 7 e contain no square in 
Bo (that is, contains only white squares and squares in Bi), and if either extremity of 70 is at a concave 
corner of U then the corresponding extremity or extremities of 7 e are at the corresponding corners of P e . 

Because of the boundary conditions on P £ , 7 E has length which is an odd multiple of e. If we remove 
7 e from P e , then what remains is a union of two disjoint polyominos Pi and P2. Let Pi be the polyomino 
which contains the base square of P £ ; then Pi is a Temperleyan polyomino. The other polyomino P2 will 
become a Temperleyan polyomino if we remove a single square s of type Bo in P2 adjacent to one of the 
endpoints of 7,, . 

Note that the union 7 e U{s} has a unique domino tiling. The number of tilings of P e equals the product 
of the number of tilings of Pi and the number of tilings of P2, divided by the probability Pr(7 E U {s}) 
that the tiling of 7 e U {s} occurs in a uniform tiling of P e . 

We can repeat this procedure on Pi and P2, cutting them apart into simpler and simpler pieces until 
we arrive at a collection of Temperleyan rectangles (a Temperleya n rect angle is an odd-by-odd rectangle 



with corners in Bo and one corner square removed). Temperley [Tem| provides us with a formula for 
the exact number of tilings of a Temperleyan rectangle (see Proposition Working by induction, 

to compute the number of tilings of P e it suffices to be able to compute the probability of finding, in 
a random tiling of P E , a tiling of 7 e U {s}. We cannot compute this probability exactly but we can 
approximate it sufficiently closely (Lemma |^). 

The region 7 e U {s} has a unique tiling which is a chain of dominos. Starting from the boundary of 
P;, let ai,a2, . . . ,<J]v be the set of consecutive dominos making up the tiling of 7,= U {s}. The en are 
laid end-to-end except for ajv which is perpendicular to the others. The probability of all the ai being 
present in a tiling of P e is a product, as j runs from 1 to TV, of the probability that a,j is present, given 
that ai, . . . , aj-i are already present: 



N 



Pr(ai, . . . , ajv) = Pr(a^ | ai, . . . , aj_i). (7) 

Suppose that ai, . . . , fflj-i are present already. These dominos form a strip running from the bound- 
ary of P e to a point in the interior of P e . The region pP' d = P E \ {a±, . . . , eij-i} is again a Temper- 
leyan polyomino, by our hypothesis that the at contain only black vertices of type Bi. So computing 
Pr(aj I 01, • . • , a,j-i) is a matter of computing the probability of a,j in a random tiling of this region P> . 

Let Uj be the region U with a slit cut out along the segment 7o([0, 2(j — l)e]), and translated by 
— 2(j — l)e so that the tip of the cut is at the origin. Then Uo — U (up to translation) and Un is a union 
of two rectilinear polygons. 

We may suppose (after applying a rotation if necessary) that the path 70 is horizontal and goes from 
left to right. For < j < N let fj be the unique conformal isomorphism sending the right half- plane 
{z : Re(z) > 0} to Uj which sends to (end of the cut), 00 to the base point bo, and has expansion 
fj(z) — z 2 + 0(z 3 ) at the origin (we do not define /o or /jv)- 

Lemma 6 Let be the function Fi (coupling function limit) on the region Uj. We have 
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where err(j) is o(max(i, j^rj)). The term in brackets on the right side is 0(e) when j is not close to 
or N . Let £ > be a small constant and K = K(e) := f/e. Then this can be written 

1 + — + em(j) + 0(e) ifj<K 
3 — 1 

1 + gS^(O) + -o(e) ifK<j<N-K, (8) 

1 + + err 2 (N - j) + 0(e) if N - K < j 

N - J 

where S^/J] is the Schwarzian derivative of and erri(x) and err^ix) are o(i). Here err\ depends 
only on whether 70 staris at an edge or at a concave vertex, and err^ depends only on whether 70 ends 
on an edge or at a vertex. The constants Co,Ci are determined as follows. If Jo starts on an edge, then 
Co — — i; if 70 starts at a corner then Co = — If Jo ends at an edge then C\ = — |; if jo ends at a 
corner then C\ — — |§. 

As we will see in the proof, the first and third expressions on the right hand side of (^) are special 
cases of the middle expression, except for the error terms. 

The probability of Lemma can be related to the change in normalized Dirichlet energy of the limiting 
average height function of Uj : 

Lemma 7 Let £ and K be as in the previous lemma. For N > j > 1 the difference in the S-normalized 
Dirichlet energy of the limiting average height function between Uj and Uj-i is Es(hj) — Es(hj-i) = 

if j < K 

if K < j < N - K, 
ifN-K<j 

where Co and Ci are defined as in the previous lemma. When j = 1 or j — N , the difference in energy 
Es(hj) — Esihj-j) is a constant depending (for j ' — 1 ) only on whether the cut starts at a corner or edge 
and (for j = N) on whether the cut ends at a corner or edge. 

These two lemmas, and a computation of the number of tilings of a Temperleyan rectangle, give us 
the main result: 

Proof of Theorem [l]. The proof is by induction on the number of cuts required to cut P E apart into 
Temperleyan rectangles. In the case P e is a Temperleyan rectangle, V = 4 and Proposition hq below 
shows that (|l|) equals the log of the number of tilings, up to an error 0(e). 

For a general P e as in the statement, let j e U {s} be a strip as explained above, cutting P e apart into 
two Temperleyan regions P' and P" . We compute the log of the probability of j t U {s} occurring in a 
tiling of P e . This probability depends on whether j t starts on an edge or at a corner, and whether it 
ends on an edge or at a corner. 

Let Oi, . . . , ajsr be the chain of dominos of j e U {s}, and let P^ = P e U {ai, . . . , a,}. By Lemma ^| we 



Pr(aj|ai, . . . ,a 3 -i) 
a/2 - 1 



48C 
irj 



+ 0(e) 



TV £ 



48Ci 
'n(N-j) 



+ 0(e) 
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have 



Pr(> U {s}) = y PrjaAa! a 3 -!) = 
(V2- 1)^ Y 1 

^log(l + ^+err 1 0) + O(e))+ £ log (l + fS^O) + i (e)) + 



j<K ' K<j<N-K 



+ log(l + ^^+err 2 (iV-j)+0( e ) 

N — K<j ^ ^ 

= E T +err3(j)+ ° (e)+ £ ^^(0) + i o ( e )+ ^ _^ +err4 (jV-j)+0(e). 

j<if •* K<j<N-K ' N—K<j 

In this expression the terms err3(j) and err4(jV — j) sum to give an error ERRo(e) log i where ERRo(e) 
is o(l). The remaining terms are equal by Lemma [7] to — 7r/48 times the change in Dirichlet energy on j e , 
up to the error terms O(e) and |o(e). However the terms O(e) sum to KO(e) — 0(£) which is negligible. 
Similarly the terms jo(e) sum to A-o(e) which tends to zero if £ tends to zero sufficiently slowly. 

Since each tile (except the last) decreases the area by 2 and increases the perimeter by 4, noting 
that — 2co + 4ci = log(\/2 — 1) which cancels the denominator of (|^), the formula (jl]) is correct up to a 
justification of the term ERRo(e). 

The error ERRo(e) is the sum of two error terms, one coming from the beginning of y e and one from 
the end. Let ERRi(e),ERR2(e) respectively be X^j<K err3 0) wnen 7« starts at an edge or a corner 
respectively. Let ERRz(e),ERR4,{e) respectively be }2 >N _ K err4(N — j) when 7 e ends at an edge or a 
corner respectively. We have ERRi(e) + ERRi{e) = ERR2(e) + ERRs(e) since when 7 e is traversed in 
either direction the same error occurs (that is, the probability of 7 6 U {s} is independent of the position of 
s as long as it is of type Bo and on the correct side of f € ; therefore 70 can be traversed in either direction, 
giving the same probability for 7,= U {s}). Moreover by Proposition [l6| ERRi(e) + ERRs{e) = 0, since 
when we cut a rectangle apart into two rectangles there is no error (or rather this error, when multiplied 
by logi, is still 0(e)). 

Therefore if 7 e begins at a corner and ends at an edge, or vice versa, the error ERRo is ERR2 — ERR\. 
When 7 E begins and ends at a corner the error ERRo is 2(ERR2 — ERRi). When we cut P e apart to make 
rectangles, each concave corner gets cut exactly once, so the total accumulated errors will be the number 
of concave corners times a fixed error ERR2{e) — ERRi(e). Setting ERR(e) — ERR2(e) — ERRi(e) 
completes the proof. □ 

As a shortcut for computing the log probability of the cut 7 E U {s}, from Lemma ^ the log of ([?]) 
is ATlog(-\/2 — 1) plus the sum over j of log ^1 + (^=eF[^ (e, 2e) — ljj , plus the error term. When e is 
small the sum over j can be replaced by an integral using the variable t = 2ej: 

_ r2e(N-l) / \ rlt 1 

Pr( 7c U {s}) = JVlogCv^ - 1) + J \-j=eF[ tI2e \t, 2e) - lj • — + ERRo{e)\og - + const + o(l). 

(10) 

Here ERRa(e) depends on whether or not 70 begins or ends at a corner, but not on U. This form will 
be useful in section |^. 



4 The probability of the next domino on a cut 
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4.1 The slit plane 

Define the slit plane SP to be the plane minus the left half of the x-axis: SP = R 2 — (— oo,0]. For 
e > define polyomino SP e to be the (infinite) polyomino obtained from P(2eT? + (0, e)) by removing 
the lattice squares centered at (—ke, 0) for all k > 0. We assume the lattice square centered at the origin 
is of type W\. The infinite polyomino SP^ is Temperleyan, with base point 6 at infinity (SP t is gotten 
from the graph H t which is obtained from 2eZ 2 + (0, e) by removing edges (— 2ke, — e)(— 2ke, e) for all 
k > 0). 

For later use, we compute the asymptotic coupling function on the slit plane SP. 
Lemma 8 The asymptotic coupling function on SP with base point at oo satisfies 

r,SP/ \ 1 1 

F + (v, z) = 



TTy/v (\/z — y/v) 

F- (v,z) 



ITy/v (yfz + ■\/v) 

Proof. The map z i— > yfz maps SP to the right half- plane RHP and oo to oo. On RHP we have 
F+ HP (v,z) = f^T) and F^ HP (v,z) = ^t^j- The result follows using the transformation rules 

(i- □ 

Let U, Uj be as in section ^. Let / = fj be as in that section. Define b = bj and c = Cj to be the 
coefficients in the expansion 

f{z) = z 2 + bz 3 + cz 4 + 0(z 5 ). (11) 

Note that b £ iR and c G R since / maps the imaginary axis to the real axis. The inverse of / has the 
expansion 

r^) = - 1/2 -^+f¥-^V 3/2 +^ 2 )- as) 



2 V 8 2 

The expansion of y/J is 



v7R = .+ ^ 2 +(f-!) 2 3 +0( 2 4 ) (13) 



and the Schwarzian derivative of yff at the origin is defined as 



a /7 rm ._ WD'" 3/(^/)"\ 2 _ 96 2 

5/7(0) - wfy~ 2 1 wJ " 3c ~~- (14) 



4.2 Proof of Lemma ^ 

Here we prove Lemma []. After a translation we may assume that the right-hand square of the domino 
Cbj-i is the square centered at the origin in el?. Then Pr(aj|ai, . . . , Oj_i) = Cj(e,2e) where Cj is the 
coupling function on Pi J ' translated as above. We must then approximate Cj(e,2e). Recall that e is a 
lattice square of type Wi and 2e is of type B\ . 

Recall from (^) that for v € Wi, the coupling function Cj (v,z) on P e satisfies AC-, (w, z) = S v + e (z) — 
6 v - e (z), for z £ H' the dual graph of H = H{P^), where A is the Laplacian on H' . Furthermore Cj(v, z) 
is zero when z is the outer face of H. In particular ACj(e,z) = $2 t (z), that is, Cj(e,z) is the discrete 
Green's function on H' centered at e. To compute Cj(e, 2e), we will use a theorem of Kesten comparing 
the values of a continuous and a discrete harmonic function near the tip of a cut: 
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Theorem 9 ([Kes]) Let D s ut be the slit disk {\z\ < 1} \ (—1,0] and let g be a continuous harmonic 
function on D s m with piecewise continuous boundary values and boundary values on the slit. Let g t be 
a discrete harmonic function on el? D D s i it with boundary values within O(e) of the boundary values of 
g, and boundary values on the slit. Then for any (k,£) G 1? , 

g c (ke, le) = X k .tg(ke, le) + a{sfe) 

where the constant Xk,e only depends on k and I. 

We cannot apply this theorem directly since in our case Cj(e, z) is not harmonic at the point z — 2e. 
However let Csp € (v,z) be the coupling function on the discrete slit plane SP £ . The function g(z) := 
Cj(e,z) — Csp{e,z) is harmonic for z £ B\ in a neighborhood of the origin in H' , including the point 
z = 2e, since the Laplacians of Cj and Csp e are equal there. Furthermore g(z) — (Dirichlet boundary 
conditions) when z is on the cut 70. In fact g(z) is discrete harmonic on all of Pj J ' except possibly where 
the interior of P e meets the negative z-axis somewhere to the left of the slit (but since Cs p e admits an 
analytic continuation around the origin, we can define g(z) so as to be harmonic on all of P e ). 

We can therefore use Kesten's theorem applied to g(z) to compute g(2e) = Ai,og(2e), where g(z) is 
the continuous harmonic function on Pg with the same boundary values as g(z). It remains to compute 
g(z) and Csp c - 

Lemma 10 On 5P e with z € Bi, and z not within 0(1) of the origin we have Csp e (e,z) = ref\ SP (e, z) + 
o(e) = TW~^ + o(e), for some constant t. 

For the proof see the appendix. Now to compute g(z), from the lemma it is the (continuous) harmonic 
with boundary values —reFf p (e, z)+o(e) on dllj. However note that the function eF^\v, z) — eFf p (v, z) 
as a function of z is continuous, harmonic, with boundary values —eFf p (v,z) on dUj. Multiplying by 
r, at v = e we must have g(z) = re(F^\e, z) — Ff p (e, z)) + o(e). 

Let Nl(0) be the neighborhood of radius L of the origin, where L is chosen small enough so that 
N L (0) C P E (j) . More precisely, if j < K or j > N - K we take L = min{je, (N — j)e} and for 
K < j < N - K take L = £ . __ 

On dN L (0), eFf p {e,z) is 0(y/J) ( see Lemma |), and eFi(e,z) is comparable to eFi P (e,z) and 

therefore also 0(^/ J). Therefore on 8Nl(Q), y^(rePi(e, z) — reFf p (e, z)) is bounded and we have by 
Kesten's theorem (for the disk of radius L) 

^g(2e) = Ai,ot J| (eF 1 (e, 2e) - ePf P (e, 2e)) + oj I), 



g(2e) = Ai j0 r [ePi(e, 2e) - eF? P {e, 2e)j + o(-). 

Now from Lemma ^ (using Pi = \(F+ — P-)) we have eFf p (e, 2e) — ^P, and plugging in the value 
of Ai.or from ( p^ ) below yields the first result. 

Let / = fj : RHP — > Uj be as in the statement. Let s(z) = f^ 1 (z). Then from the transformation 
rules (^), (^|) and (g) we have 

7T{S(Z) — S(v)) 

P_ [V, Z) - 



n(s(z) + s(v)) 

Consequently 



ePi(e, 2e) = - 



S '(e) , s'(e) 



ty \ s(2e) - s(e) s (2 e ) + s (e) 
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Using the expansion ( |l^ ) and (|l4|), and the facts that b £ iR, c € R this reduces to 

TV 07T 

So then _ 

Pr(a J ia 1 ,...,a j _i) = y/2 - 1 + Ai, r^S , - v /^(0)e + o(|). 

Plugging in for Ai,qt gives the center result in (^). 

When j is close to we must replace e in Theorem |^ by j since j is the combinatorial distance 
(number of lattice points in H) to the boundary of the region. We have 

(a J |ai ! ...,a 3 -i)=v / 2-l + A 1 , r^5 A /^(0)e + o(i). (15) 



Pr 

67T J 

A similar expression holds when j is close to N. 

Now when j is small or close to N the Schwarzian derivative of fj is blowing up in a standard way 
which is independent of U up to an error O(l): 

Lemma 11 When j is small the germ of fj at the origin is independent ofl] (depending only on whether 
or not the cut starts at a corner or on an edge) up to an error 0(1). Similarly when j is near N the 
germ of fj at the origin is independent ofU up to an error 0(1), only depending on whether the cut ends 
at an edqe or a corner. 



The proof is in the appendix. Below when we write O z (l) we mean an analytic function in z each of 
whose coefficients is O(l) (as j — » or j — » N). In particular z 3 O z (l) is an analytic function whose 2-jet 
vanishes. 

Recall that fj is normaliz ed as in (|ll| ). If 70 starts at an edge of U, then when j is small fj is 
approximated by the map 2^/ ejz 2 + (ej) 2 — 2ej which maps RHP to {x + iy \ x > —2ej} — [— 2ej, 0]. 
That is, from Lemma [ll] we have 



fj(z) = 2y/ejz* + (ej) 2 - 2ej + z 3 O z (l) = z 2 - — + z 3 O z (l). 
Therefore Sy/fc(0) = -jfj + O(l). Plugging this into (|l|) yields 



Pr(a j \a 1 ,...,aj- 1 )=V2-l - ^ o 7 1 +o(-) + 0(e) 

where the o(j) term is independent of U: the term o(i) only depends on the error in Kesten's theorem 
as applied to the coupling function on {x + iy \ x > — 2ej} — [— 2ej, 0]. 

If 70 starts at a corner of U, then fj is approximated by a map of the form 

f(z) = C + C'(z + 2iB)Vz~iB + z 3 O z (l), 

for constants C, C' and B (this follows from the Schwarz-Christoffel formula: the derivative of / is a 
constant times ^*_ iB f° r some real B). The conditions that f(iB) — — 2ej and f(z) — z 2 + 0(z 3 ) 
determine C,C' ,B and a short computation gives 

fj(z) = z --z-j=z - -T-. +0(z + Z Oz{l), 

whence 



16 



A similar argument holds when j is near to N. If 70 ends at a straight edge of U, then near the origin 



fj is approximated by the map 



c" 



+ z 3 O z (l), 



Vz 2 + B 2 

which maps RHP to the region shown in Figure |^a. Again C, C' , B are determined by /(oo) = 2e(7V — j) 
and f(z) = z 2 + 0(z 3 ), and a calculation gives S^/J](0) = + 0(1). 

Lastly, if 70 ends at a corner of U, then near the origin fj is approximated by a map of the form 



fj(z)=C- 



c 



Vz - iB(z + 2iB) ' 

which maps RHP to the region shown in Figure |§|d. Using /(z) = z 2 + 0(z 3 ) and /(oo) = 2e(jV — j) 
determines C,C' ,B to give 



/(*) 



- 1 

z + 



2iV3 



9^J 



4cj 



7 z 4 + 0(z 5 ) + z 3 O z (l), 



which gives S^f~(Q) = - Jfj + O(l). 







,2(N-j)e 







2(N-j) £ 



Figure 3: When the cut ends at an edge (a) and at a corner (b). 

This completes the proof. □ 

To determine Ai ; oT, we use Proposition |l^ below. When we cut an m x n rectangle R into two 
rectangles Ri and R2 with a horizontal cut of length m, the sum of the logs of number of tilings 
of Ri and R2, minus the log of the number of tilings of R, which is the probability of the cut, is 
m log(\/2 — 1) — i log m + O(l). On the other hand from the above proof this is 



ni ) 
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This gives 



\l,OT= -j= . (16) 



5 The change in the Dirichlet energy 

We prove here Lemma [?]. The computation is in two steps. We first show that the change in the <5- 
normalized Dirichlet energy only depends (up to an error o(e)) on the 4-jet of / at the origin, that is, 
only depends on b and c of (|ll|). Then we do an explicit computation of the Dirichlet energy for a 
particular 2-parameter family of regions, whose corresponding maps / have 4-jets covering every possible 
value of (b, c) . 

5.1 Dependence on 4-jet 

Since in section |2.4| we computed the dependence on 8 of Eg, we can assume without loss of generality 
in this section that 8 is small compared to e. 

Let hj denote the limiting average height function for Pj and hj+i be the limiting average height 
function for Pj+i. 

Lemma 12 Up to an additive constant we have hjofj — -^/rralog f'j(z). In particular hj = hjofjofT 1 — 
-f/mlogt/- 1 )'^). 

The proof of the first statement follows from the definition of hj in the last paragraph of section |2.S| . 
The second statement follows trivially from the first. 

Let gj be a harmonic conjugate on Uj of hj, so that hj := hj + igj is analytic on Uj. Similarly let 
gj+i be a harmonic conjugate for hj+i on Uj+i. 

Since the boundary of Uj is a subset of the boundary of Uj+i, the change in normalized Dirichlet 
energy can be written (recall that Uj is the region Uj minus a (^-neighborhood of its vertices) 



hj+i dgj+i - f hj dg 3 = * (h j+1 dg j+ i - hj dgj) + f hj+i dg ]+1 (17) 
au' , Jau' Jau' Jx 

(hj+i - hj)dg J+ i - * hjd{g j+1 - gj ) + <f> h j+1 dg j+1 (18) 
au'. Jau'. Jx 

3 3 

where X — dUj +1 — dU'j is the path which consisting of the four pieces: dNg(0) U [—8, 2e + 8] UdNs(2e) U 
[— 8, 2e — 8]* where the superscript * is a reminder that the second segment is traced in the reverse 
direction. 

We'll show that each of the three integrals of (j^|) depends (up to controlled errors) only on the 4-jets 
of fj and f J+1 . 

The first integral in (hj) can be estimated as follows. On the bound ary o f Uj we have hj+i = hj, and 



for each corner c of Uj+i we have dgj+i = O(S) on dNs(c) (see section 2.4). It remains to consider the 
boundary dNg(0), where is the tip of the cut in Uj (which is not a corner of Uj+i). The a;-axis divides 
dNg(G) into two parts. Since gj+i is smooth on each half, when 8 — > the integral of (hj+i — hj)dgj+i 
on each half tends to zero. In conclusion the first integral in ( |l8| ) tends to zero with 8. 

To compute the second integral in (|ri|), we show that we can replace the path of integration by a 
smaller path of radius O(e); on this new path we will be able to replace dgj and dgj+\ with their 4-jets. 

Now hj+i — hj is the harmonic function which (after choosing an appropriate base point) is zero on 
dUj, and 1 and —1 on the upper and lower side of [0, 2e], respectively. On dUj we therefore have 

9i+i - 9i = Im(/i J+ i - hj) = -i(hj+i - hj) 
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since the real part of this function vanishes. 

The second integral in ( p^ ) can then be written (using hj = — ^lmlog(f^~ 1 )' (z) from Lemma |l2[) 

f hjd(gj+i-g j )=--Imcf log^YMH) (h j+1 - hj) dz, (19) 

where the ' refers to the ^-derivative. This integral is the imaginary part of a contour integral. The 
integrand is analytic on U' — [0, 2e]; in particular it is analytic on U' minus a neighborhood of the origin 
of radius 3e. Therefore we can replace the path of integration by a path which winds around the boundary 
of the slit disk D sHt = {z e U' :_\z\ < 3e}. 
Now on D a ut we have (see (J12|)) 

log(/r 1 )'(^) = log(^= - \ + ~/^+zO z (l)) = log((/r { 1 4} )')(2) + * 3/2 O z (l) 

where fjj 4 \ is the 4-jet of fj at the origin (7 is a constant depending on b, c only). Similarly 

logC^+r)'^) = Iog( 2 ^_ 2e - ^ + 7J+i>/^2i + 0(z - 2e)) 
= log((/7 + 1 1 , {4} )')(^) + {z- 2 e ) 3 / 2 O z _ 2e (l), 

where recall that since we are using coordinates in Uj, fj+i(z) = 2e + z 2 + frj+iz 3 + Cj+iz 4 + . . . . 
These equations give, on dD 8 u t , 

dg 3 = -^Re (dloga/r^)') + v^O,(l)d2) = dg JtW + ^feO z {l)dz (20) 

and similarly dgj+i = d<?j + i,{4} + y/eO z (l)dz. 

Having found dgj and df/j+i, the integral ( |l9| ) becomes 



hj{dg j+ i - dgj) = * (^,{4} + z 3/2 2 (l)) (d^ +lj{4} - dg jt{4} + y/eO(l)dz) 

JdD slit K ' 



(21) 



Since hj is bounded, ^flO(X) § hj^ 4 ydz = 0(e 3 ^ 2 ) and similarly on the path of integration the only 
singularity of dgj+x,{4\ or dgjj 4 \ occurs near the origin, where the functions are at worst 0(z ) so 

j z 3/2 O z (l) (dg j+1 , {4} - dg J>{4} ) = j z 1/2 O z (l)dz = 0(e i/2 ). 

Thus we can pull the errors in (^]) out of the integral, giving 0(e 3//2 ). Therefore the integral ( |l9|) only 
depends on bj, Cj, bj+i, Cj+i up to an error Q(e 3//2 ). 

Lastly we estimate the third integral in ( |lq ) . As was the case for the first integral, on the two halves 
of Ns(0) the integral of dgj+\ tends to zero with 5. Now on Ns(2e) we have dgj+i = O(S), and hj+\ is 
constant (1 or —1) on the two "sides" of [0, 2e]; so 

r2e-S rS 

hj+idgj+i = o(l) + / dgj +1 - j dg~ +1 (22) 

X JS J2e-S 

= o(l) + g J+ i(S) + - ffi+i(2e - 5)+ - g j+1 (2e -5)-+ g j+1 (5y 

where the + and — refer the the two limits from above and below the axis, and the o(l) tends to zero 
with 5. Using gj+i(z) — ■|Relog(/ :i+ 1 1 )'(z), we see that the values at 5 and 2e — 5 of gj+i depend only 
on the 4-jet up to error 0(e 3//2 ). 
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We have shown that the change in energy depends only on the 4-jet of fj and fj+i- Now when j is 
not within if of or N, 6j+i = 6 3 + O(e) and c 7 + i = Cj +0(e) by Lemma ftU Changing bj+i and Cj+i by 
O(e) changes the second and third integrals and by at most 0(e 3 ' 2 ); so when K < j < N — K 
the change in energy in fact depends only on the 4-jet of fj . 

When j < K or N — K < j, the energy depends on the 4-jet of both fj and fj+i- However in these 
cases both fj and fj+i are independent of U up to 0(1); in the next section we'll see that the energy 
depends only on j (through the Schwarzian derivative of \J fj) up to error o(l/j). This completes the 
proof that the change in Dirichlet energy only depends on the 4-jet of f at the origin, up to terms in 
o(e). 



5.2 Computation for a specific family of functions 

At this point we could simply work out the integrals ( [jl] ) and (^ij), but it is less computationally painful 
to do an explicit calculation for a particular family of functions. This will also give a more accurate 
estimate of the energy when j is near or TV. 
For p, q £ iM. and {0,p, q} distinct, define 

f p , q (z) = 2 x p- f uj^—ldu. (23) 

V p Jo V u - 1 

If — ip > > — iq and jpj > |g| this is a map from the RHP injectiv ely o nto the region U PtQ shown 



in gray in Figure |4|a; this follows from the Schwarz-Christoffel formula [Ahl| (/ maps the RHP to the 
polygonal region with angle 2n at the origin, 3n/2 at f(p), n/2 at f(q), and 2n at oo. If — ip > > — iq 
and \p\ < \q\ this is still a well defined map which is not injective on RHP (although it is locally injective). 
This lack of injectivity is of no consequence to us. When p and q are on the same side of the origin the 
image is as shown in Figure ^3. Again the map may or may not be injective according to whether or not 
p is between q and 0. 

The 4-jet of f p , q at the origin is 

2 fl l\2 3 /3 1 1 \ Z 4 



For any / of the form ([n]) there is a p, q for which f P:q has the same 4-jet at the origin as /, provided 
that b / 0: it suffices to take p = ^ziJU 2 an< ^ 1 ~ i6c+9i^ • Note that when p, q take these values, two 
of the points {0,p, q} are equal only when 6 = (the cases when one of p, q is infinite are allowed). The 
case 6 = will be dealt with later. 

The integral in (B3|) can be explicitly evaluated, giving 



/*,«(*) = iV~ (V(^ - PW - - 2P + 69) - 2(p 2 + 2pq - 3q 2 ) log (V^P+ V^q)) - C 2 



where C2 is chosen so that f P , q {0) = 

- /2 

4 V V 



° 2 = l\f^(VM(- 2 P + 6<Z) - (P 2 + 2pq - 3q 2 ) log {-(y/p - y/q) 1 



We have explicitly 



f(p) = -]q(2p -6q) + ljl(pl + 2pq~ 3q 2 ) log ^ + ^ 



4 4V pV J \</q-</p 



and 



/(«) = -i«(2p - 6.) + i/f (P 2 + 2 M - ) log 
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f(p) 



f(q) 



f(p) 



f(q) 



Figure 4: The image of the RHP under the function j VA is shown in gray. 



Thus 

/(p) -/(?) = x\/f (p2 + 2M " 3a2) (25) 

where we chose the sign of the square root and branch of log to correspond to the situation of Figure Q 
Suppose — ip > > — iq. Then the 5-normalized Dirichlet energy for the region U p , q = f PA (RHP) is 

§jj, hdg where ft is on the vertical boundary from — ioo to f(q), — 1 on the horizontal boundary from 

f(q) to 0, 1 on the boundary from to f(p), and 2 on the vertical boundary from f(p) to oo (note that 

h has different values on the two "sides" of [f(q), 0]). 

To compute Eg, we pull h back to RHP. The preimage of Ng(f(p)) is to first order a disk around p 

(that is, it converges to a round disk when 8 — » 0); let 5 P denote its radius. Similarly let S q and <5o be 

the radii of the preimages of the disks around f(q) and /(0), respectively. Let Soo be defined by: l/<5oc 

is the radius of the preimage of the disk of radius 1/5 around 0. 

Then Eg is (recalling that g is constant up to lower order terms on the boundary of the ^-neighborhood 

of the singularities) 

r-i& rp-i&p pt/8oa 

Eg = - d{g°f P ,q)+ %°/p l9 ) + 2 / d{gof p , q ) 

•J q + i5 q J 5q Jp+ig p 

= -2g(f M (8o)) + g{f P , q {q + iS q )) - g(f P , q (p + i5 p )) + 2g(f p , q (l/S 00 )). (26) 
From ([m]) we have 8q — S up to higher order terms. For z close to p we have from (^) 

/*,«(*) = UAp) + - pf' 2 + o{(z - pf) 

3 vp- q 
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and so up to higher order terms we have 

2/3 



and similarly near q 



U,M = f 9 ,M + 4^qVq~p(z - qf /2 + 0(z - q) 



*.= (7i/!r^)- (28) 



,4V q iVl-P 

For large z we have fp >q (z) = 2. / |z + O(l) which implies f P ,q(z) = J ^z 2 + O(z) and 



S °° = ( S ][~) 1/2 - W 
Plugging these in to ( [2d] ) with g o f PA = — ^ Re log f^ q gives the Dirichlet energy to be 

+ >'-(^l(^ S )-^(^(l)" 4 ) '»» 

= -log^ + -logf 8 ^ , 8 . (31) 



7T \ \j q / J Ti I I /P 1 

\ 4 V 9 1-Jq-V , 



3tt b 3tt b V4 8 g 19 |p- g| 

By symmetry we get the same energy when — ip < < — iq. A similar calculation shows that the 
same energy is obtained in the remaining cases when p and q are on the same side of the origin. 

When we extend the cut by a small amount 2e, the new region has a new uniformizing function 
of the same form / p ', 9 ', where p and q' are defined by changing p and q so that f(q) and f(q) each 
change by — 2e. Let dp and dq denote the changes in p and q for an infinitesimal e. We must first have 
= d(f(p) — f(q)). This relates the change in p to the change in q. From (Eq) we have the equation 

= d(f(p)-f(q)) 



r(2p + 2«) + (p"+2pg-30(s^) ) d P+ ( a/z(2p- 6 9) + — (P + 2 W-3g") ) d 9 , 



«n=-2fy +:jw+3< cw 

p \P + 6pg — 15g 2 / 
Since . / -(p 2 + 2pq + 3q 2 ) does not change, we have d(f(p)) — 



-ldp+(Sq4p)d< 1 +lAp+^-^) ( ( A + ( "F^F + 1 ^ ^ 



l&pq 2 



(p 2 + 6pq — 15q' 



-dp. 
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Now the change in E$ when d(f(p)) 



= -2c and d(/(p)) = d(f(q)) is 



-2e 



dEs dp 
dp d(f(p)) 



4(5p + 7q) (p - g) (p 2 + 6pq - 15g 2 



np(p 2 + 6pq — 15g 2 ) 
_ e(5p + 7g)(p-g) 



-16pq 2 



2%p 2 q 2 



Finally the Schwarzian derivative of \J f P . q is 



3 
4g 2 



1 

2p~q 



1 

4p2 



1(1 
4 l g 



P 



(5p + 7g)(p-g) 
16p 2 q 2 



which is j times the change in E$ . This completes the proof when j £ [K, N — K] except in the case 
6 = 0. For the case 6 = see the function f q of section |5.3.1 below. The case j £ [K, N — K] is dealt 
with in the next section. □ 



5.3 Beginning and ending of a cut 

As above, the change in energy near the beginning of a cut only depends on the germ of / near the tip of 
the cut. Near the beginning of the cut, the germ of / only depends on whether or not the cut starts at 
a corner or at an edge (Lemma [n]) . We computed the limiting functions / when j is near the beginning 
or ending of the cut in the proof of Lemma 

5.3.1 beginning on an edge 

Suppose first that the cut begins on a straight edge of U. We compute the change in E$ between the 
time when there is no cut to the time when the cut has length 2ej. 

Let h q be the function on RHP — [0, q] whose boundary values are on the y axis, 1 on the upper 
boundary of the cut [0, q] and —1 on the lower boundary of the cut. Up to an additive constant, this is 
the height function on U near the beginning of the cut. 

We can compute the <5-normalized Dirichlet energy of h q in the same way as we did the functions f Ptq 
of the previous section. Specifically, the map from RHP to U q = RHP — [0, q] is 



/,(*) = ^2qz 2 + q 2 



q 



y/ + OAz% 



and the Dirichlet energy is 



2fl°/ 9 (5o)-3°/«(J-|+5i)-flo/ 9 (-^-| + 5 2 ), (32) 

where (to first order) f q maps the (^-neighborhood of the origin to the 5-neighborhood of the origin, and 
the neighborhoods of ±\/ — f °f radius 8i to the 5-neighborhood of —q. 



As before So = S 1 ^ 2 ; one computes Si = 82 = — —, Plugging into (p2|) with go f„ = — -Relog f'„ 



gives 



^ = -logi + -logg+-log2. (33) 
When q — 2ej, i.e. U q — RHP — [0, 2ej], this gives the energy change due to the first j steps of the 

8e , 



cut. The change in Dirichlet energy between j and j + 1 is ^ which is —^S^/fjiO) as required. 
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Note that if we set 5 = e in (|33|), the e- normalized energy due to the beginning of the cut is 

— logj + const + 0(e). 

TV 

In particular when j = K = £/e this is — log 1 + O(l). 
5.3.2 starting at a concave corner 

In case the cut starts at a corner of U (of angle 3tt/2), the change in energy can be computed in a similar 
manner. Again we compute the change in energy between the time when there is no cut to the time 
when the cut has length 2ej. 

Let TQP be the three-quarters plane, TQP = {(x, y)\x > or y > 0} and define U q = TQP — [0, q]. 
We first compute Es on TQP. We can take the average height function h on TQP to be 1 on the 
negative x-axis and on the negative y-axis, that is, h — const + ^lmlog(z). Then E$(h) on TQP is 
§ hdg — g(S) — <?(|) and up to higher order terms g(8) = -2- log ~ and g(j) = — j- log ~. So the energy 

The map from RHP to U q is 



/,(*) = 2V~q ^ -^L = 2V~q qf 2 + 2q( W - qM 

Jo V w -i V 3 / 



At the origin we have 

Now Eg = § u: hdg — —go fq{j—) — 3 ° fq(q + iS q ) + 2g o f q (8o). Again <5q = S 1 ^ 2 and near oo, 



x-i ft 1 \ 4 v /3 q 3/2 
s = AT-) = — 5 — <*°o ■ 

Near g we have f(z) = /(g) + Aq^J—q^Jz — q + 0(z — q), so 

|W=9l<^ /2 = 5. 

Plugging all this in we have 

Es ( hq) = I log ( 2V~ q ( ^ ) + J log - 2 1 log 



4^5 

14, 1 20, 
= it log t + it l°g 9 + const. 

The difference in energy between U q and TQP is then 

10 , 1 20 , 

ir- log 7 + — log g + const. 

Now /(g) = which is — 2ej when g = y^-- Plugging this value of g in gives the energy due to the 

IT 1°S 7 + T" lo s( e j) + const - ( 34 ) 
When j changes by 1 the energy changes by -^r, which is again — — S'-y/^'fO). 

Setting 8 — e in (|34j) gives the contribution to E t from the beginning of the cut as log j + const + 
O(e). 
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5.4 Ending on the interior of an edge 

A similar computation holds in this case, with f q as in the proof of Lemma ^. However there is a shorter 
method, obtained as follows. 

The change in Eg near the end of a cut only depends on the local structure of the average height 
function near the end of the cut. Therefore the energy can be obtained from the known value of the 
energy for a rectangle. Indeed, cutting a rectangle into two rectangles with a single edge, the change in 
E e due to the beginning of the cut, plus the change due to the ending of the cut, plus the contribution 
from the "central" terms, must equal total energy change which is 21 log | + c onst + Q(e) (see fell)). The 



contribution at the beginning of the cut is ^ log i + const + 0(t) (see section 5.3.1 ), so the contribution 
at the end of the cut is 

18 , 1 ~, s 48/3, 1 , 
— log — h const + U(e) = — — log — h const + U[e 

TV £ 7T \ 8 e 



The constant is independent of U by Lemma [H]. 



5.5 Ending at a concave corner 

The energy for the end of a cut ending at a corner of U can be computed from the previous case. Given 
an ell-shaped region, cut it into two rectangles with a single cut beginning on an edge and ending at 
the corner. We can compute the change in Dirichlet energy by starting the cut at the concave corner 
and ending at the edge; the contribution from the start of the cut starting at the concave corner is (cf. 
(0)) !vr 1°S e const + 0(e) and the contribution from the end of this cut is log I + const + 0(e) 
by the previous paragraph. The sum of these must equal the contribution for starting from the edge 
(^ log ^ +const + 0(e)) and ending at the concave corner. Therefore the contribution due to a cut ending 
at a concave corner is 

46, 1 . 
— log — h const + U(e) 
3n e 

This completes the proof of Lemma m. 
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= ™ ( lo s 7 + const + °« 



6 The case of rectangles 



Let P be a (2m — 1) x (2n — 1) rectangle with a unit square removed at the lower left corner. By 
Temperley's trick [ rcrr| , domino tilings of P are in bijection with spanning trees of the m x n grid graph, 
rooted at the lower left corner. 

The number of spanning trees of the m x n grid, which by Kirchhoff [ Kir | is the product of the 
non-zero eigenvalues of the Laplacian, is 



l-2ros| — 



2cos[ ^ 



n 

where j (resp . fc) runs from to m - 1, resp. to n — 1. The log of this formula was asymptotically 
computed in | DD |: 

Proposition 13 (|DD|) The log of the number of spanning trees of an n x m rectangle is 

^52 + (m + n) log(V2 - 1) - \ log(m) + \og{n{e~ 2 ™ /m )) - - log(2) + 0(— ). 
7r Z 4 mn 
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Here G is Catalan's constant and r\ is the Dedekind eta-function 



1/24 



n a-A 



The area of P is (2n - l)(2m 
formula may be rewritten 



G 

7T 

where 



Area(P) + 



G_ 

2tt 



io g (V2-i; 



1 = 4nm — 2n — 2m, and the perimeter is An + Am — 4, so this 



Perim(P) - i log(m) + \og(n(e- 27Tn/m )) + C + 0(m" 2 ), (35) 



C = log 



2 

2 5/4 

T+71 



+ 



2G 



6.1 Dirichlet energy for a rectangle 

Let U be the rectangle [0, |] x [0, J] in C, with base point at the lower left corner. The average height 
function h for U is, up to an additive constant, the harmonic function which is 0, 1, 2, 3 on the lower, 
right, upper and left boundaries respectively. 

The funct ion h has an explicit expression in terms of the the Weierstrass elliptic function p(z) — 
pi,ir{z) (see [Ahl : this is the doubly-periodic function with periods 1 and ir and a double pole at each 
point of the lattice). 

Lemma 14 Up to an additive constant we have 



Proof. From [Ahl we have 



h(z) = 
p'( z ) 



-Imlog p'(z) 



E 



(z — w) 3 ' 



where F = Z + irZ. Note that p'(z) is real when z G [0, |] or z £ [0, ^] + i5 and pure imaginary when 
z £ i[0, |-] or z € i[0, ^] + i. Furthermore on a fundamental domain for T, p'(z) is zero or infinite only 
at the corners of U (see [ Ahl | ) . So the argument of p'(z) is constant on each edge of the rectangle U. At 

; ) and so — Imlog p'(z) changes by +1 at 

□ 



a corner Zj, Zj 7^ 0, we have p'(z) — Cj(z ~ Zj) + 0((z — 
each corner of U when going counterclockwise around U. 



Let U' be U minus the 5-neighborhood of the four corners of the rectangle. Then 



E s (h) = / hdg = 



1/2-6 



0dg + 



(1-Kt)/2-»5 



ldg + 



1/2+iS 



r/2 + S 



(l+i-r)/2-S 



■ + 



ir/2-iS 



3dg + 0(S) 



where the contribution on the boundaries of the neighborhoods of the corners is O(S). Since when 5 is 
small g is essentially constant on the neighborhoods of radius 5 of the corners, this energy is 

E s (h) = - S (i - 5) - g{\ + il - 5) - g(i T - + S) + 3g(S) + 0(6). 



At the origin, p'(z) = — % + 0(1) and so 



log p'{Se l 



•log 



p e 3i 



+ const + 0(6) 



2G 



°Icg±-™+*Icg(-2) + 0(0. 

TT TT TT 



Near the other three corners of U, 



2 2 

- log p'(<5e' e -Zj) = \og{c 3 5e ) + 0(6) 

TT TT 

where Cj = p"(zj). Using g(z) = ^ Re log p'(z), the energy is 

2 18 1 

E s (h) = (log(ci<5) + log(c 2 5) + log(c 3 <5)) + — log(-) 

tv no 

= — log- log C1C2C3), 

7T TT 

where ci, C2, C3 are the derivatives p"(zj) at the three other corners of U. 
Lemma 15 We have 

I/O \12 / -27tt\24 

C1C2C3 = 2^ 27V > ) • 

Proof. From the differential equation 

(p'( z )f = Hp(?) - ei)(p(z) ~ e 2)(p(z) - e 3 ) 
(where ej = we obtain (upon differentiating the logarithms of both sides) 

2 p"(z) _ p'(z) | p'(z) | p'(z) 



y(z) p(z)-ei p(z)-e 2 p{z)-ez 
from which we get 

p"(zi) = 2(ei - e 2 )(ei - e 3 ), 
with similar expressions for C2 , C3 . Their product is 

C1C2C3 = 8(ei — e2) 2 (ei — e3) 2 (e2 — e-s) 2 



and by [Apl this equals 



The energy is therefore 



A 1 / \12 / -27TT\24 

-A=-(2tt) 77(e ) 



.log- log -(2tt) T)(e 7TT ) 

7T 7T \ 2 



□ 



This is the 5-energy for a | x § rectangle. 

Let (7 a , /3 be the rectangle [0, a] x [0, where r = /3/a. We can scale U a> /3 by ^- to get {/. As a 
consequence the above expression is the 2a<5-energy for U a ,p. To get the 5-energy for U a ,p, substitute 
5 /(2a) for 8 in the above expression. We get 

E s (U a>0 ) = HI log H£ _ 2 log C 1 (2ff )«„( c -*-)« > ) . (36) 
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Proposition 16 Let U be a rectangle; let P t be a sequence of Temperleyan rectangles of area A t and 
perimeter Perim e approximating U . Then the log of the number of domino tilings of P e is 

^ + ( f + 1 -^4^-) ^ - U) + C + { e% 

ne 2 \2n 4 J e 48 

where r^ie, U) is the e-normalized Dirichlet energy of the average height function (given by putting e for 
8 in ) and C is a universal constant. 

Proof. If U is an a x (3 rectangle, this follows from (^) upon setting 2m — 1 = a/e, 2n — 1 = j3/e. 

□ 



This proposition tells us why the 'natural' choice of S in Eg(h) is S = e or some constant multiple of 
e. In fact there is a multiple which makes the constant C in the above formula vanish, although we will 
not need this accuracy. 



7 Loop-erased random walk 



The loop-erased random walk in a finite graph G has the following simple description. Let bo, &i be two 
vertices of G. Take a simple random walk starting from 6o and stopping at b\ . Erase from the path its 
loops, in chronological order. That is, if there is a loop, erase the first loop (the first time the path comes 
back to the same vertex twice). For the new path, if there is still a loop, erase the first loop, and so on. 
The remaining path is a simple path from bo to &i. 



There is a well-known connection between the LERW and spanning trees [Pern]: in a uniformly chosen 
spanning tree on a region P C Z 2 , the unique arc (branch) from a to b has the same distribution as the 
LERW from a to 6 on the same region. Pemantle also showed that on an infinite domain such as 1? 
or th e upper half of Z 2 , the LERW is well-defined (despite recurrence of the simple random walk. See 



7.1 LERW and dominos 

Let P be a 2n — 1 x 2n — 1 square Temperleyan polyomino with base square bo on the center of the 
right edge. Recall how Temperley's trick works: a domino tiling of P corresponds to a spanning tree of 
the graph H(P) which is rooted at bo in the following way. Each domino covering a vertex v of H(P) 
has white vertex covering the center of an adjacent edge e of H{P). In the associated spanning tree the 
outgoing edge from v points in the direction e. The union of these edges forms a spanning tree with all 
edges oriented towards the root at &o- 

Let Q be the region obtained from P by removing a single black square b € -Bo and a single white 
square w. 

Lemma 17 When b is on the boundary of P, domino tilings of Q are in bisection with spanning trees of 
P for which the branch from b to bo contains the edge w (traversed in either direction). 

Proof. Let T be a domino tiling of Q. Let b' and b" be the neighboring squares of w which are in 
H . Temperley's trick assigns to a vertex x G H an outgoing edge e if and only if xe is a domino of T. 
Temperley's trick, applied to T, gives a set of directed edges of the graph H such that each vertex has 
exactly one outgoing edge, except for bo and b which have no outgoing edges. Furthermore the union of 
these edges is a forest, i.e. each component is a tree. There are exactly two components to this graph 
since each component is rooted at exactly one of bo or b. 

We claim that b' and b" are in different components of this forest. To see this, first note that the 
directed branch from b' cannot pass through 6", for otherwise the path from b' to b" followed by the 
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segment from b" to w to b' would be a lattice path in H enclosing an odd number of squares of Q |^ and 
so could not arise from a domino tiling of Q. Similarly the paths from 6' and b" cannot end both at 60 
or both at 6 because their union would contain a closed lattice path from b' to 6". 

So the paths from b' and b" end one at bo and one at b. On the path which ends at 6, shift the 
dominos along it by one square towards b, and add an extra domino from w to its freed neighbor. This 
makes a domino tiling of P in which the tree branch from b passes through w (shifting dominos by 1 
along a tree branch has the effect of changing the direction of the edges on that branch). 

This process is reversible: from a tiling of P whose branch from 6 to bo passes through w, shift the 
dominos by one up to w to get a tiling of Q. □ 



When b is not on the boundary of P, it is possible that one of the paths from b' or b" winds around 
b, or that these two paths both lead to bo, enclosing b. So the lemma does not hold in that case. 

By Lemma to prove Theorem ^| it suffices to be able to count the number of tilings of Q, or rather, 
to compute the ratio of the number of tilings of Q to the number of tilings of P. 



7.2 Region with a white hole 

First we compute the coupling function on Q. Since Q has a hole (the hole to) which does not enclose 
the same number of black and white squares mo dulo 2, the weighted adjacency matrix A of Q' (the dual 



graph of Q) is not a Kasteleyn matrix for Q. In [Kas2|, Kasteleyn describes how to redefine the weights 



of A to get a Kasteleyn matrix in this case (rather, he describes how to get weights for general planar 
graphs). One way to get a Kast eleyn matrix for Q' is the following. Start with the weighted adjacency 



matrix A = A(Q') of section 2.2. Take a path 7 of vertices in Q from the outer boundary to a vertex of 
w (7 is a path of faces of Q 1 ). Every edge in 7 crosses an edge V1V2 of Q' . For each edge in 7, change 
the sign on the corresponding matrix entries A V1V2 and A V2V1 . The matrix A' with these new signs is 
a Kasteleyn matrix for Q' . In particular it has the property that its d eterm i nant i s the square of the 



number of matchings, and its inverse is the coupling function Cq for Q [Kas2, Kenl]. 

Let Q' be the double cover of the graph Q' , branched around the face containing w. That is, Q' is 
defined by the property that a closed path in Q' lifts to a closed path in Q' if and only if it winds an 
even number of times around the face w. Similarly we define the double cover H of H as follows. Recall 
that w is an edge of H. Remove w from H and let H be the double cover of H — {w} branched over 
the face which contained the edge w. The double cover of H' , the dual of H , is similarly defined (and 
denoted H'). 

For each fixed v, the coupling function Cq(v, z) lifts to a discrete analytic function Cq on Q': let z' 
be a lift of z and v' a lift of v; then define Cq(v' , z') to be ±Cq(w, z), with the + sign when v' , z' are on 
the same sheet of the cover (in the sense that there is a path from v' to z' which does not cross a lift of 
7) and — sign otherwise. 

For a point x £ Q' let x' and x" denote its two lifts. 

Let v € Wi and z € B\. For lifts v' of v and z' of z, the real part of Cq(v',z') is harmonic on H' 
except when z — v±e or z = w±e (as in (Q), with extra singularities at w ± e). So letting G denote the 
Green's function on H' , and v± = u±e, w± = w±e, it must be that KeCg(v' , z') is a linear combination 
of the Green's functions G(v'±, z'), G(t>^, z'), G(w'±, z'), G(w'±, z'). By the antisymmetry under changing 
sheets and (^), ReCc^w', z') is of the form 

G(v' + , z') -G(vL , z') - G(v'i , z') + G(v'l ,z')+ ai (g(vj'+, z) - G(wL , 2')) + «2 (<S« , z') - G{w'L ,z')), 

for some (real) constants ai,«2 which depend on v. These constants are determined by the following 
condition on Cq: the harmonic conjugate of ReCcj must be zero at the four points 60,61,60,6". These 



1 Prom the Euler formula for a disk F — E + V = 1, for any simple closed path in 7? the sum of the number of faces, edges 
and vertices strictly enclosed is odd. 
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four conditions give only 2 linear equations for cti, ai since the value at b' ' is by construction the negative 
of the value at b' and similarly for 6". 

We need to show that the coefficients ol\ , a-z converge as e — > 0, and compute their limit. To show this 
requires two steps. First, we will show that the Green's functions G(y, z) converge to the corresponding 
continuous Green's functions, and moreover the terms \(G{w' 1 ,z)-G{w'{,z)) and ±{G(w' 2 ,z)-G(w'l,z)) 
converge to the derivatives of the continuous Green's functions. These derivatives have simple poles at 
the origin. Then we will show that there is a unique pair of analytic functions Fa, Fi with the required 
properties having simple poles at the origin. 

The convergence of the Green's functions on H' is standard: As e — > 0, if y does not converg e to 
then G{y , z) has the form G(y' , z) = § log \y' - z'\ + | log ± + const + 0(e/\z' - y'\) (see e.g. [3pJ]). 
This is enough to conclude that as e — > 0, \{G(y — e, z) — G(y + e, z)) tends to a harmonic function with 
zero boundary values and a singularity at y of the form Re^^^y . This is the derivative -§^g{y, z) of the 
continuous Green's function g(y,z). 

From the symmetry G(y,z) — G(z,y) we can conclude the same when y tends to as e — > 0. In 
conclusion ±(G(w[, z)-G(w" , z)) and \{G(w' 2 ,z)-G{ui2,z)) converge to harmonic functions with simple 
poles at the origin (and single- valued harmonic conjugates). 

To show that cti and oti converge, we claim that it suffices to show unicity of the limit of Cq . That is, 
since for each e, ai and «2 are solutions of a linear system whose coefficients converge (being functions of 
the Green's function derivatives), either in the limit the linear equation becomes singular (in which case 
there is non-unicity of the limit) or it does not, in which case the solution is unique and the solutions for 
finite e converge to this solution. 



7.3 Unicity of the limit 

Rather than work on a square region, it suffices to work on the unit disk with the white hole at the 
origin: the transformation rules (Q, (|B|) allow us to move the result back to the square Q. 

Let U be the unit disk with marked points bo, 6i on the boundary. Let U be the double cover of U 
branched over the origin (the map f(z) = z 2 maps U to U). The lifts to U of the asymptotic coupling 
functions Fo(v, z), Fi(v, z) on U have zeros at z = ±\fba~, ±\/b\, poles at z = ±«, and simple poles at 
z = 0. (A priori the poles at the origin may have residue 0, in which case they are not poles at all). 
These functions also are antisymmetric: for v, z £ U we have Fo(v, z) — —Fo(v,—z) = Fo(—v,—z) and 
the same for Fi. Furthermore, since Fq and F\ are respectively pure imaginary and real on the boundary 
of U , they extend by Schwarz reflection to meromorphic functions on the entire Riemann sphere, with 
additional simple poles at ±l/w and oo (the reflections of ±v and 0). In particular they are rational 
functions. Since we know exactly the location of their (six) poles, and four of their zeros, it remains to 
find the other two zeros. The antisymmetry under z i— > — z and symmetry under z i— > 1/z implies that 
the two other zeros are on the unit circle: otherwise the orbit of the zeros under these two symmetries 
would have four elements. Thus the functions have the form: 



F (v,z) = c 



Fi(v,z) = ici .,.,„.,_.,. 

{z 2 - — ir)(2 ■* — jr) 

for real co,ci and where 63 = 63(11), (respectively 64 = 64(11)) must have absolute value 1. Now 63 and 64 
are chosen so that the residues of Fo and Fi are real at z = v. The constants co = co(v),ci = ci(v) are 
real and chosen so that the residue of each function at z — v is —. 



Vbo \ f z _ y6i \ f z_ 63 

60 Z ) \^/bl z J \b 3 z 



(z 2 — v 2 )(z 2 — V 2 ) 

; Vb~o~ \ ( _z Vbi \ f z_ _ 64 

60 z J V^6i z J V64 z 
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We can explicitly solve for 63 as follows. The residue of _Fb at z = v is 



Co - 



( v _ VM.\ I v _ V^i \ ( j> 

V v'Eo J \ Vbl v / \ b 3 « 



2v(v 2 — v 2 ) 



This is supposed to be real; equating this and its complex conjugate yields 



( V j ( V Vbl \ f V 



yv^7 ^ J \ " J \ b 3 

2v(v~ 2 — v 2 ) 
Solving for 63 gives b% = Xq /Xq , where 



iTST) l 6 ^- ssr 



2u(u 2 — v 2 ) 



X = 



bo 

bo v 



v 



— )v + «V6o - —7= 
« / V vVbo 



vs/bi 



Vy/bi 



Similarly one finds b\ — —X1/X1, where 



60 

fb~o v 

A lengthy calculation yields 

F+(v,z) = 



V y/bi\ f r- 1 

7=-— + [vVh-^-= 
lb\ v J \ v^/bo 



vvb 



Vyfbl 



■iz 



z v&i 
Jb\ z 



7T Z — V 



(which is meromorphic in both v and z as expected) and 



4z 



F-(v,z) = 



« vfci 



/bi z 



7r(l — z 2 w 2 ) I vy/b~o — 



«V&1 — 



/biv 



(anti-meromorphic in v as expected). 

The map from U to U is /(z) = \fz. Using the transformation rules (^),(|5|) we have on the original 
region U 



F-(v,z) 



2{z - b ){z - b x ) H[ 
tt(z — v)(v — bo)(v — bi) V z ' 

2{z - b )(z - bi) 



7r(l — 2w)(?;6o — l)(«6i — 1) V z 

The fact that the transformation rules apply in this case follows from the fact that the results F± 
have all the required properties of the coupling function limits (and are the unique functions with these 
properties) . 



7.4 On the RHP 

The map f(z) = maps the RHP to the unit disk, sending 1 to 0, to —1 and 00 to 1. Let RHP° be 
the RHP with a 'white hole' at 1. The limiting coupling functions on RHP° with zeros at 65 := / _1 (&o) 
and be, := are 

p rhp°, n_ 2(2- b 5 )(z- be) 
f+ l u > z ) 



tt(z — v)(v — bs)(v — be) V z 2 — 1 ' 
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F *«r°( v z) = 2(z-b 5 )(z-be) v 2 
1 ' ' it(z + v)(v + b 5 )(v + b 6 ) \ z 2 
We will assume ba = oo ; in this case the formulas on RHP° become 



j^RHP" , s 2(z - b 5 ) ' 

n(z — v)(v — OS) V z — 1 



F MF °(v) = - ? 2(z : t5 ' 1 / T ^. (38) 

V ^ 7r(z + + 65) V Z 2 - 1 V ^ 

7.5 The number of tilings of Q 

Let U be the K x K square U = [0, K] x [— if/2] C RHP. Suppose Q — Q e approximates U as 
e — *• 0, and suppose that the holes b,w,bo of Q 6 , converge to points of the same name b, w,bo G 17. 

Since we are concerned primarily with the growth rate of the LERW, to simplify the calculations we 
will make the further assumption that b,w £ U are within distance 1 of the origin, bo is on the right edge 
of U and K is large. Using Lemma [ll] we can then approximate to error 0(1/ K) the coupling function 
limits F U , FY for points near the origin by the (much simpler) coupling function limits on RHP — {b, w}. 

Suppose that w is on the z-axis at ^-coordinate a > and let b = f3i where j3 is real, |/3| < 1. 

As in Lemma ^ let P be the region Q with the holes b, w filled in. Let S C P be the chain of 
horizontal dominos from the origin to w. 

Notationally, for a region A let N(A) denote the number of tilings of A. The ratio of the number of 
tilings of Q to the number of tilings of P is a product of three terms 

N(Q) = N(Q) N(Q - S) N(P-S) 

N(P) N(Q-S) N(P-S) N(P) ' [ ' 

each of which we can now approximate. 

7.5.1 The first ratio 

The inverse of the first term in (pE|), N(Q — S)/N(Q), is the probability of S occurring in a tiling of Q. 
This is computed in a similar fashion to the proof of Theorem |IJ Let ai . . . , a m be the set of dominos in 
the chain S. Then the probability of cij occurring given 01, . . . ,ctj-i is given by (see Lemma ^|) 

where Uj is the translate of the region RHP — [0, 2ej] (with a white hole at w and black hole at b), 
translated by — 2ej so that the tip of the cut (2ej) is at the origin; Uj corresponds to the polyomino 
Qj = Q — {a\, ■ ■ ■ ,a,_i}. The limiting coupling function F Q Uj on Uj can be computed from (^) and 
(p8|) . Indeed, let t = 2ej € [0, a]; then the map 



' Va 2 - t 2 

maps the region Uj to RHP and sends the singularities a — t to 1 and i/3 — t to 65 = f(i/3 — t). 
A computation using @ and (||) gives 

eF^(e,2e) = ^ (f(e), f(2e)) + F^° (f(e), /(&))) 

V2(5t 2 - a 2 )e 3/2 
7T 47rt(a 2 -i 2 ) 1 ; ' 
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As a consequence 

Pr{aj\ai,... ,Oj) _ (5t 2 - a 2 ) 

y/2-1 4t(a 2 ~t 2 ) + ° [e) - 

When we sum the logs of Pr{a,j\a\, . . . ,%-i) for j running from 1 to m we get (see equation (|Io|)) 

log Pr(S) = m \og(V2 - 1) + jH"* |p^y ■ | + o(log e) + o(l) 
= mlog(\/2 - 1) + - log - + i log a + o(l) + o(loge) 

O € O 

where the o(log e) term is independent of a and (3. 

Surprisingly, there is no dependence on f3 to first order. 



7.5.2 The second ratio 

Since 6 and w are on the boundary of the polyomino P—S, the second ratio in (^) is exactly \C p ~ s (w, b)\, 
where C p ~ s is the coupling function on P — S. 

The coupling function on P — S can be computed using (^), Q and (^) and the map f(z) — 
v/(5 + a) 2 — a 2 which maps RHP — [0, a] (translated by —a so that the tip of the cut is at the origin) 
to RHP, and oo to oo. We must compute eF p ~ s (e, —a + i/3). We have 

eF ( r S (e, -a + i/J) = £/'(e) ( /( _ q + _ /(e) + /( _ Q + ^ + /(e) ) = ^=f + 0(e)- 



7.5.3 The third ratio 

The computation of the third ratio in j39| ) is similar to the first ratio, except that now there are no 
singularities at b,w. We may use the coupling function on RHP, given by (^|), to compute the coupling 
function on RHP — [0, t] (again shifted so that the tip of the cut is at the origin). The appropriate 
function is f(z) — y/Jz + t) 2 — t 2 and we arrive at 



The log of the probability of S occurring is then 

mlog(V2 — 



e dt 



= mlog(V2- 1) - ~ log a- \ log - + o(loge) +o(l). 
o o e 



7.5.4 The product 

Now combining these three results, the log of the ratio ( ps| ) is 

log(-R) = -mlog(v / 2-l)-i log i-i loga+log — ^|EL=+m log(^-l)-i log a-i log -+o(loge)+o(l) 
8 e 8 7r\/Q 2 + /3 2 ° 8 e 

3 11 1 
= -7 log - + - logo - -log(a 2 + 1 ) +o(loge) +o(l), 

where the o(loge) is independent of a,/? (the o(l) depends on a,/3). 
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This is the log of the ratio of the number of tilings of Q and the number of tilings of P. For finite K 
there is an additional additive error term 0(1/K). 

Letting a + i(5 — re' e in polar coordinates, the probability that w is on the LERW from b to 60 in P e 

is 

/e \ 3/4(1+0(1)) . 

(I) cos(0) 1/4 (l + O (l)) 

where the o(l) in the exponent is independent of 6. 

The expected number of edges on the LERW from 6 to bo and which are at distance < R from b is 
then the integral of ^£^ 3/ ' 4 ( 1+0 ( 1 " cos(9) 1/ ' 4 , times the number of edges per unit area times the area 
form r drdO, as r runs from to R. This gives (.S) 5 ' 4 ( 14 ~ W_ This completes the proof of Theorem^. 

8 Open problems 

There are a number of places where our arguments could stand improving, or where there are interesting 
avenues for further research. We list some of the outstanding ones here. 

1. There remains in Theorem ^ the error term ERR(e) which is somewhat annoying. It seems rea- 
sonable to suspect that this term is in fact O(l) but we don't have any method of computing this 
term at present. 

2. Is there an extension of Theorem [j] to the non-simply connected case, as well as to the case of 
surfaces of higher genus? A more general version of Lemma ^| as well as its proof should be easy 
given the coupling function _Fb. We don't know at present an easy extension of Lemma to this 
case though. 

3. Can one find a more natural proof of Lemma [j]? There may be a proof using invariance properties 
of the Schwarzian derivative but we did not see it. 

4. Theorem |l| could be generalized to regions with polygonal boundaries which are horizontal, vertical, 
or have slope ±1: this is because there is an exact formula for the determinant of the Laplacian on 



a triangular region {(x,y) : x > 0, y > 0, x + y < n}, see [KPW 



More generally, one can ask about the relationship between the term — 75-02( 6, U) of (|l|) and the 



(^"-function regularization of the determinant of the Laplacian. It appears from OPS | that they are 
the same up to a multiplicative constant. Is there a simple explanation for this fact? 

Is 5/4 the almost sure growth exponent for the LERW? The major open question about the LERW 
is: Is there a scaling limit of LERW, and if so is it conformally invariant? 



9 Appendix 

9.1 Green's function for a slit plane 

We prove here Lemma [To| 

Proof. As noted in the first two paragraphs of the proof of Lemma |^, we must compute the 
discrete Green's function G(0, 2) on 1? — (—00, —1], that is, the function of 2 harmonic on 1? — (—00, 0] 
with Dirichlet boundary values on the boundary (— 00,— 1], which satisfies AG(0, 2) = So(z) for z G 
1? — (— 00,— 1], and asymptotically G(0, z) — > when |zj — *■ 00. 

Let fn(z) be the harmonic function on Z 2 — (—00, —1] whose boundary values are for z G [— n, —1] 
and 1 for 2 G (— 00, — n — 1]. Let g n (z) = f n {z) — f n +i(z — 1). By definition, g n (z) is harmonic on 
Z 2 — (— 00, 0] with boundary value / n (0) at 2 = and boundary value for 2 G (— 00, — 1], Therefore g„ 
is a constant times the desired Green's function G(0, z). This constant is 1 over the Laplacian of g n at 
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0, that is, — l/(/n+i(0) + /n+i(i) + fn+i(— i)) (note that f n {z) is harmonic at and f n +i(z — 1) is zero 
at z — and z = — 1). So we have 

G(0,*)- Uz)-f n+1 {z-l) 



/n+l(0)+A+l(*) +/«+!(-*)' 

The asymptotic values of /nfz) for large n and 2 are ilmlog gy Theorem lol applied to the 

71 1 — iWz/n L - ' 

function f n (z), we have /«(0) = Ai.o-^= and / n (i) = fn(—i) = Ao,i Also when »i 3> |z| 3> 1 we 
have ^ 

- /n+i(z - 1) = — -/= + 0( ). 
Therefore when \z\ 3> 1, and in the limit as n — » oo, we have 

f 



G( °' 2) Uai,o+2v / 2A ,J v^ + ^ )_ 
Scaling everything by e (replacing z with z/e) completes the proof. □ 



9.2 Local Riemann mappings 

We prove Lemma [Tl| This is essentially a weaker version of a lemma of Hayman [Hay, Lemma 6.6]. 

First consider the case when j is small and the slit starts on an edge of U. For q > let g q be the map 
g q {z) = \J 2qz 2 + q 2 — q. Then g q maps the RHP injectively onto the translate by —q of RHP — [0, q}. 
When q = 2ej, g q is the standard limiting form for the fj \ we must show that the derivatives of fj differ 
from those of g-ztj by O(l). There is a constant C (depending on U but independent of j) such that for 
each j sufficiently small g^jfj maps a neighborhood of the origin in RHP to the half-disk Bc(0) CiRHP. 
Now g^jfj extends by Schwarz reflection to an injective map from a ball around the origin into Bc(0); 
since the derivative at the origin of g^jfj is 1 this ball has radius comparable to C. Since g^jfj is 
bounded on this ball, the Schwarz Lemma implies that its derivatives at the origin are all bounded as 
well. We conclude that g^Jjfj = z + z 2 O z (l). 

Now fj = g2tj(z + 2 2 O z (l)) and since <?2ej is independent of U, the germ of fj at the origin is 
independent of U up to O(l). 

The same argument with a different g q works in the case where the cut starts at a corner. The 
argument at the end of the cut requires a simple modification. □ 
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